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We present analytical expressions for the current-current correlation function in graphene for 
arbitrary frequency, wave vector, doping, and band gap induced by a mass term. In the static limit 
we analyze the Landau (orbital) and Pauli magnetization, as well as the Lindhard correction which 
describes Friedel and Ruderman-Kittel-Kasuya-Yosida oscillations. In the nonrelativistic limit we 
compare our results with the situation of the usual two-dimensional electron gas (2DEG). We find 
that the orbital magnetic susceptibility (OMS) in gapped graphene is smeared out on an energy scale 
given by the inverse mass. The nonrelativistic limit of the plasmon dispersion and the Lindhard 
function reproduces the results of the 2DEG. The same conclusion is true for the Pauli part of the 
susceptibility. The peculiar band structure of gapped graphene leads to pseudospin paramagnetism 
and thus to a special form of the OMS. 

PACS numbers: 



I. INTRODUCTION 

Graphene, which was first isolated in 2004,1^ is 
the name of a monolayer of carbon atoms that are 
arranged in a hexagonal lattice. It differs compared to 
most two-dimensional systems by its relativistic energy- 
momentum relation and its nontrivial spinor structure,!^ 
originating from the two-atomic Wigner-Seitz cell, and 
has remarkable electronic properties.- As a conseque nce, 
various effects like the anomalous quantum Hall effeclP^ 
or the Klein tunneling have been discovered. In 
1956, the orbital magnetization of two-dimensional 
graphite had already been calculated,'^ indicating a 
strong diamagnetism in the undoped hexagonal lattice 
which was confirmed by recent experiments.'^ 

In the present work, we study the response of the 
system to an electromagnetic potential in terms of the 
current-current correlation function. Similar studies have 
already been performed recently regarding the density- 
dcnsity response of masslesa^Eo] and massful-'J'' Dirac 
fermions, the Ha ll con ductivity in the presence of spin- 
orbit interactionJi^Iini including self-energy and vertex 
corrections,'^ and current-current correlations in the 
absence of a mass term. ^^'^^ Here we will generalize those 
results to the case of massive quasiparticles by taking into 
account a mass term which breaks sublattice symmetry, 
leading to a gap between the valence and the conduction 
band. It can occur due to diffe rent mechanisms, 
including intrinsic spin-orbit couplin^^^E^ (with a gap 
of graphene placed on a suitable substrate 

{Eg = 0.26eV)^ or adsorption of molecules (with a 
gap of several electron volts) ^ The current correlator is 
related to the polarization function, which was discussed 
earlier.'^f^ Its limiting behavior determines the orbital 
and the Pauli magnetization, the plasmon spectra, and 
the screening of electric or magnetic impurities. Without 
the mass term, the Landau magnetization is infinite for 
intrinsic graphene (i.e. zero chemical potential, fi — 0) 
and zero for extrinsic graphene (/x ^ 0) P | i5 | i6 | 22 | 2 3] 
the Pauli part vanishes for the former and is finite 



for the latter case. As gapped graphene is similar to 
the two-dimensional electron gas (2DEG), we investigate 
in the nonrelativistic limit, i.e., the limit of single- 
particle energies just above the band-gap parameter, and 
compare our results to that of the 2DEG.'^ 

The particular features of the density-density and 
current-current correlation functions in graphene 
compared to the standard 2DEG rely on the coupling 
of the orbital degrees of freedom to the sublattice 
of pseudospin. On the other hand, semiconductor 
systems involving coupling to other internal degrees 
of freedom such as the physical electron spin have 
also been analyzed recently with similar aspects. As 
examples, we mention studies of the dielectric function 
of semiconductor 2DE Gs with various types of spin-orbit 
coupling terms,'^^'2ll two-dimensional semiconductor 
hole systemSjl^^ and p-doped bulk semiconductors.'^ 
Moreover, the dielectric function of graphene taking into 
account the full honeycomb lattice structure (but not a 
mass term) was analyzed recently in Ref:^^. Analytical 
expressions for the polarizability of graphene with finite 
width of Landau levels, temperature, and mass term can 
be found in Rei^. 



This paper is organized as follows. After introducing in 
Sec.|ll]the model Hamiltonian and pertaining quantities, 
we present in Sec. |III| analytical expressions for the 
longitudinal and transversal current-current correlation 
function. 



In Sec. IV 



we focus on the static limit 
and determine the orbital and Pauli magnetization. 
Moreover, we include many-body effects via random- 
phase approximation (RPA). In Sec. |vj we study the 
effect of an increasing mass term on typical quantities 
like the magnetic susceptibility, Friedel oscillations, and 
the plasmon spectra, and we compare the results to 
the 2DEG. We close with conclusions in Sec. IVII In 
Appendix [A] one can find details of the calculation of the 
transversal susceptibility, while Appendix [B] comments 
on the relation between current and density response. 
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II. THE MODEL 

The atoms in graphene are arranged in a honeycomb 
lattice, where each unit cell contains two carbon atoms. 
The effective Hamiltonian near the corners of the 
Brillouin zone K/K', including a mass gap as well as finite 
doping, is given by, using standard notation 



T^x i my /i — ravp 



(1) 



where we have introduced an electromagnetic vector 
potential via tt = hk + A. The upper (lower) sign 
refers to the point K (K'). The field operator is defined 

'ay 

operators of the Bloch states in the two sublattices. 
Concentrating on the K point [upper sign in Eq. ([T])], 
the eigenvalues and eigenspinors at zero vector potential 
A = Q are given by 



by *fe = 



where Ofc and bk are the destruction 



E± {k) ^ fi± \J {hvpk)'^ + {mvj^ f, 



lx±(fe)> 



1 ± 



V2 U./lT 



y/{hk)^+(mVFf 



with k = \/ + k'y. The current operator follows from 



SHo 
SA 



(2) 



where a are the Pauli matrices. Equation ^ is, up 
to vp, equal to the pseudospin operator. The electric 
current can be connected with the vector potential via 
the correlation function Xj^j^ i defined by the Kubo 
produclPl] 



A{t),B{0) 



Our system is rotationally invariant and the current is 
thus a linear combination of a purely longitudinal and a 
transversal part {q = \q\): 



where f{E) is the Fermi function, and g counts orbital 
and spin degeneracies ((7 = 4 in graphene). 

The orbital magnetic susceptibility is given by the 
static transversal part of Xjj ■ 



Xorb 



e2 x^- 
— hm 



(9,0) 



(4) 



Because of the continuity equation, idtPq — q ■ jq, 
the response to a scalar potential, i.e., the polarization 
function, is included in the current-current susceptibility. 
In graphene, this leads to the following relation (see 
Appendix iBlPEl 



Xppiq,^) = q^Xjjiq,^) 



hA 







(5) 



The second term on the right-hand side was calculated 
in Ref.^^ and reads 



1 

hA 



QJq,P-q 



where D is a cutoff parameter, which is usually chosen 
to be of the order of the inverse lattice constant.!^ Note 
that the commutator is independent of the mass. 

For the following, it is essential to distinguish the cases 
mv"^ > p and p > mvp. In the first, intrinsic case, the 
Fermi energy lies between the two bands, while in the 
second, extrinsic case, the Fermi energy lies either in the 
conduction or in the valence band. From here, we will 
omit the spatial indices and use the notation Xjj = Xj^j^ ■ 



III. RESULTS 

We restrict our discussions without loss of generality 
to positive frequencies w, chemical potentials /i, and 
All other cases follow from x!jf^ {q,—uj) — 



mass TO. 

x^!^i<i,^) 



and by observing that the results only 



depend on the absolute value of /i and to^ 



-Xiiil,^) 



+ U, 



For a nonintcracting system, Xj^j^iQi^) is given by 

f{ExAk))-f{Ex.Ak + q)) 



A 



^ huj- 

Ai ,X2,k 



Ex,{k) - Ex,{k - 

X (XAi {k)\a^\xx^ (k + q)) 
X {x\2 (k + q)\a-f,\xx, (k)) , 



A. Intrinsic case 

In the intrinsic case, only transitions from the valence 
into the conduction band contribute. As described in 
the last section, the longitudinal part, i.e., q — qx, can 
be obtained from ^ and the density response given in 

Ref.^^. The longitudinal part was also directly calculated 

i) + iO by the authors in order to check relation ^ for finite 
TO. Because of the similarity to the transversal case, we 
restrict details of the calculation, given in the appendix, 
(3) to the latter. The results are: 
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FIG. 1: The different regions related to the 
imaginary part of the current-current correlation 
function for the extrinsic case of mvp/i.i = 0.9. See 
Table |I] for the definitions of 1A-5B. 



TABLE I: Definition of the different 
regions in the g-cj-plane related to the 
imaginary part of the current-current 
correlation function in the extrinsic 
case, cf. Eq.(|9l). 



lA: 
IB: 



2A: 

2B: 

3A: 
3B: 
4A: 



4B: 
5B: 



VV"jy \-2 . V 

I < 2kF A yJihvpqY + 4(mu|,)2 < huj 



< 



±/i =F \/{hvF)'^{q - kpY + {mvl)^ < huj 
fj. + ^/{hvFy{q - kFY + {rnvj-Y < Tu^ 



<li+ ^{hvFYjq + kFY + {mvlY 
fiLo < -ti + ^{hvFyjq - kFY + {mv%y 



nuj>n+ y/{hvFY(q + kFY + {mvlY 
< hvFq 



q > 2kF A ^/{hvFqY + 4:{mvj,)^ < fiw 
< M + ^{hvFYiq - kFY + {mvlY 



hvFq < huj < (hvFq)'^ + 4(mu|.)2 



h? (a;2 - [vpqf) 



(ihwf - {hvpqf ~ (2m«|)') , (6) 



Re 



g {D - 2mvl,) 



gmvpq^ 



5W 



6 [vpq — Lo) arccos 



2mvp 



6 {lo — Vpq) arctanh 



2'mvp 



Re 



9 {vpq — w) arccos 



2mvp 



(w2 - (vpqy) J 
6 {lo ~ Vpq) arctanh 



2'mvp 



h^uP- - (vpqf 



(7) 



(8) 



where d ix) denotes the Heaviside step function. transitions contribute (see above), while the second takes 

into account intraband transitions. Like in the intrinsic 
case, the longitudinal part is related to the density- 
B. Extrinsic case density susceptibility via ([s]): 

We have two contributions for the extrinsic case. The 
first one is the undoped part where only interband 
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Here we used the shorthand notation 



Gn 



x\/ Xq — ± (2 — Xq) arccos 



X\ X 



with xo = y 1 



Cq ± (2 — Xq) arccosh 
Cq ± (2 — Xq) arcsinh 



X 
Xq 
X 

xo 

X 



nmvpqy-uj-^) ^ and the regions (lA)- 



(5B), defined in Table |l] and Ref.". The chemical 
potential is defined as fi — \/ {hvphpY + [mvp)'^. The 



above functions are one of the main results of this work. 
In the absence of a gap, we recover previous results !^^^ 
Figure [T] illustrates the structure of the regions related to 
the imaginary part for the specific choice mvp/^ = 0.9. 



IV. STATIC LIMIT AND MAGNETIC 
SUSCEPTIBILITY 



In the static limit, the purely real transversal 
susceptibility is given by 



T,ext 



(9,0) = 





gmv% 








gvpq 




8T:h 



hvpq y 



/ 2mvF \ \ 



(2mvF 



1 - 



2kF 

q 



f 2fi \ 

arccos 

\hvFqxoJ 



(11) 

^[q~2kF). (12) 



while the longitudinal part vanishes, except for the 
constant term in front. Figure [2] shows the function 



-^n^(g,o) = ^(xJj(g,o) 



gP ^ 

S,TTff 



5 



for different values of a = mvp/^. 

We now insert the above functions into Q. 
intrinsic part, 



Aoro 



56 



127r c^m 



The 



(13) 



is finite and diamagnetic. Compared to the gapless case, 
the orbital magnetic susceptibility (OMS) is smeared 
out on a scale of 1/m. This broadening of Xorb also 
occurs in thepresencc of disorderjl^^l as well as for finite 
temperature!^ From (12), one can see that n-^'^^* {q, 0) = 



for q <2kF and thus 



Aoro 







which is the same as for ungapped graphcnc. The same 
result, namely, 



Xorb 



12'KC?m 



(14) 



was obtained earlier by energy considerationsP^EH The 
limit m = reproduces previous resultsPE^^^ 



Xorb 



9 2 

ge^Vp 



5(m) 



The expressions for the magnetization given above are 
only valid for the noninteracting system. A simple way 
to include many-body effects is via the random-phase 



approximation!^ The OMS in RPA is given by 



^n^(9,o) 



56 



9^0 1 



-nT(9,0) 12TTch 



[rnVf 



where eg is the background dielectric constant. One can 
see that screening effects do not change the Landau part 
of the magnetization. Without a mass gap, the RPA 
result 



A.oro 



Xorb 



yields to a renormalization, but the OMS remains infinite 
and zero, respectively. The situation changes, however, if 
one includes interaction effects in first-order perturbation 
theory beyond RPA, leading to paramagnetic behavior in 
doped graphene sheets.!^ 

The spin correlation function of a noninteracting 
system equals the density-density susceptibility.'^il The 
Pauli contribution to the magnetization follows from the 
limit 



XP = ^J■B i™xs,s,(g,o) , 



where fiB — 2m c Bohr magneton and mo is the 

electrons bare mass. The static polarization read^^l 



Snhypq 



XTn (9,0)^ 



5M 



1 - 



'1- 



arccos 



2kp 



C2mvp 
^(2mvFY + 



[hvpqY — (2mv 



2hvpiJ. q 



arccos 



2^ 



'{hvpqY 



{2mvp)'^ 



(15) 

'{q~2kF) 

(16) 



The Pauli part vanishes in the intrinsic case, refiecting 
the absence of states on the Fermi surface, while the 
extrinsic part is finite: 



V. NONRELATIVISTIC LIMIT 



A. Magnetic susceptibility 



XP 



56^ 



STTm'^c'^Vp 



mvp^ 



(17) 



Figure [2] displays the static polarization for different 
ratios a = mvp/^. The limit a — > 1 reflects the 
nonrelativistic case. 



The static transversal correlation function for the 
2DEG,l2il 



T,2DEG 



feO) = - 



91 
247rm 



1- 1 



9{q-2kF) 



leads to the OMS 



~2DEG 
Xorb 



247rmc^ ' 
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FIG. 2: (Color online) Static current correlation function for (a) extrinsic [Eq. (12)] and (b) intrinsic [Eg. (Ill] 
graphene for different ratios a = mvp/fi in units of 2i^~- Static polarization for (c) extrinsic [Eq. (16)] and (d) 



intrinsic [Eq. ( 15 )] graphene in units of ^^^a^a 



where 5 is a degeneracy factor. As described in the last 
section, the Pauli contribution to the total magnetization 
is given by the static polarization functiori^^ 



2DEG 
Xpp 



(9,0) 



l-9{q^2kF)\ l- 



2kF 



and leads to 



~2DEG _ ,,2 .9"^ 
Xp — i^B 



Figure [3] displays the function 



Xt. 



ge^m 



_^ T,2DEG/ 

^2q2 ^jj 



56 



127rmc-^ 



l--9{q-2kF 



1 - 




for the special case m — rriQ. Its limit g — )■ determines 
the total magnetic susceptibility: 

„2 



~2DEG 
Xtot 



ge^m 



3 V m / / 12^0^171 



(19) 



Expanding the graphene Hamiltonian ([!]) in the 
limit n/mvp <^ 1, and eliminating the lower spinor 
component, one finds^ 



Ha = 



2m 



:g*^J■BB 



(20) 



where o* = 2^ is the effective Lande factor, k is 
dependent on the valley, i.e., k — 1 for the K point 
and K = +1 for K'. Equation (20) is the well-known 



Hamiltonian of the 2DEG, including a Zeeman term 
which changes its sign by interchanging the two valleys. 
This Zeeman term, however, has nothing to do with the 
splitting of the energy levels due to the real spin, but 
is a truly band structure effect. Because of this, the 
second part of ( 20 ) is denoted as the pseudospin Zeeman 
term.E21 If we neglect states with negative energies, then 
the susceptibility associated with Ha is that of (18 1, 
while the magnetization is given by (19). At the same 



time, the OMS of extrinsic graphene, i.e., for ji > mvp, 
including only intraband contributions, is given by the 
paramagnetic term 



-intraband 
Xorh 



127rc^m ' 



(21) 



which means that the OMS of gapped graphene without 
hole states reproduces the total susceptibility of the 
2DEG, i.e., the sum of the Pauli and the Landau part. 
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FIG. 3: (Color online) (a) Intraband part of the transversal current correlation function in graphene for different 
ratios a = mvp/ ^. (b) Sum of Pauli and Landau contribution in the 2DEG for the special case m = toq. Both 

quantities are given in units of ]^2^c^m ■ 



Additionally, (17) describes the Pauli part due to the real 

^2^2 

spin. In the nonrelativistic limit /i ~ mvp H — gnf"' ^'i- 
(17) reads 



Xp 



ge^ra 



(22) 



which is just the result of the 2DEG. Note that (22) is 



true for extrinsic graphene with and without interband 
contributions. 



B. Priedel oscillations and plasmon dispersion 

Because of the divergent first derivative of the 
Lindhard correction (i.e., the static polarization) at q — 
2kp [see Fig. [2](c) for a > 0], Friedel oscillations in 
gapped graphene behave differently compared to the 
gapless case, where the first derivative is finite but the 
second diverges [see Fig. ^c) for a — 0]. The system's 
reaction to charged impurities is described bjff^ 



'^'totaiir) = - 



sin (2fci?r) 



and is similar to the induced spin density 6m{r), which 
describes the interaction between magnetic moments, 
e.g., due to magnetic impurities: 



dm (r ) oc — 



{2nr 

g [mvF^y ■ {2kF)^ sin(2fcFf) 
{2ttvf)'^ fifi {2kFr)^ 

fc2 2 

cases, the nonrelativistic limit reproduces the result of 



fc2 2 

Here, an = 2^ is an effective Bohr radius. In both 



the 2DEG,1^ 



Q 



sin (2fcj?r) 
(2kFrf 



dm{r) cx 



gm{2kF) sin(2/ci?r) 
{2Trfh {2kFrf ' 



while the massless case yields a different power laWj^^l 

sin {2kFr) 



^totai{r)/Sm{r) oc 



{kFry 



The long wavelength limit of the longitudinal 
susceptibility determines the dispersion of the collective 
modes. '^^ While plasmons are absent in intrinsic 
graphene, their dispersion for the extrinsic case read^^ 



Hujp (q) 



ge^^i 



\ 2eo 



1 - 



mvF^ 



In the nonrelativistic limit, this can be approximated as 



hujp(q) 



I 2tt (ehf n 



meo 



■ q 



(23) 



which equals the 2DEG result^* and particularly shows 
the same -^/n density dependence in contrast to the n^/^ 
behavior of m = 0!^ 



C. Behavior near the threshold u — vpq 

The longitudinal current correlation function for 
graphene without bandgap is singular at a; = vpq- 
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This solely results from the linear dispersion relation. 
In gapped graphene, however, the singularity vanishes 
and the response quantities discussed in this work are 
smeared out on a scale of 1/m. This is in accordance 
with the Lindhard function of the 2DEG,^' which is not 
singular at the threshold cj = vpq- However, both the 
imaginary and the real part of x^f^'~^ ^-re finite, while 
Im {xpp(g, w_f9)} vanishes in graphene [see Eq. (j9| for 
4A and 5B] and is thus in contrast to the 2DEG result. 
Furthermore, the real part at uj — vpq also differs from 
the result of the 2DEG. 



case. As gapped graphene is formally quite similar to 
the 2DEG, we studied the nonrelativistic limit of the 
magnetization, the Friedel oscillations and the plasmon 
dispersion. We have demonstrated that all of these 
quantities, which follow directly from the transversal or 
longitudinal current correlation function, can reproduce 
the corresponding 2DEG results (e.g., the n^/^ density 
dependence of the plasmon spectra or the l/r^ decay law 
of the Friedel oscillations), but with one particularity, 
namely, the pseudospin Zeeman coupling. 



VI. CONCLUSIONS 
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Appendix A: Details of the calculation of the 
transversal susceptibility 

In this section, we present details of the calculation of 
the transversal part of the current-current susceptibility. 
At zero temperature, Eq. ([3| can be written as 



Xw(g, ^) = + - with 



"47r2 



(rnvF)^ -h^k^ {1-2 sin^ipk) -h'^qk cos {ipk + Vq)\ ,.2 , 2^2 it, i,\2\ 



E{k)-E{k + q 



1 



1 



huj T E{k + q) + E{k) + iO huj ± E{k + q) - E{k) + iO 



The plus (minus) sign corresponds to Ai = A2 (Ai = — A2). (pk is the angle between k and the x axis. For the 
longitudinal case {q = qx), we obtain cos {tpk + •Pq) = costpk, whereas the overlap for the transversal part (q = qy) 
is given by cos {ipk + 7r/2) = — smipk- We now set for brevity h = vp = I. 



We define the expression 



1. Imaginary part 



A 

8nJ 





d^k 



e 1 



fc2 (1 — 2 sin^ ip) + qk sin (p 
E{k)-E{k + q) 

^2 - ^2 _ 2TUjE{k) 



2qk 



G 



6(tuj- E{k) + aE{k + q] 



for Lu > q, 



Gt. ( ^Vfc^+m^ ^l foi.q>a;. 
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where the functions < and Xq are defined in Sec. Ill 



The imaginary part for intrinsic graphene is given By 



Im 



X 1 



4m^ 



,2 _ q2 



16 



In the doped case, the upper integration hmit is not a cutofi^ parameter but is the Fermi wave vector, and we thus 
need a distinction of cases as to whether k is cut off by kp or not. For this, we define different regions, which are 
given in Table |l] and displayed in Fig. [l]in Sec. |III[ The intraband contribution to the imaginary part reads 



a,T=±l 



IGtt 



9\/\^^ - g^l 



lA 
IB 
2A 
2B 
3A 
3B 
4A 
4B 
5B. 



The addition of the intrinsic part yields the final result given by Eq. ^. 



2. Real part 

The easiest way to find the real part of the intrinsic susceptibility is by using the Kramers-Kronig relation: 



D 



X 



X Im 



9{q — (jj) arccos 



2m 



{io — q) arctanh 



2 m 



Note the cutoff-dependent part on the right-hand side. The real part of the extrinsic system is given as follows: 



Re 



g \ - / ,9. 1 



T = ±l 



,2 _ 1.2 



fc^ (l — 2 sin^ ip^ + qk sin ip 
E{k) ■ E{k + q) 



1 



1 



Lj + E{k) -TE{k + q) uj - E{k) + TE{k + q) 

k2 



2iTq^ 167r ^ ^^^^ 



2uj 



-aE{k) < 



+ / 2E(k)+cruj 



Gi 



+ ( 2E{k)+cru} 



G, 



.+ / 2E(k)+<ju \ 
V 9 J 



for q > u! 

for uj"^ > Am? + q^ 



for q^ < (jj'' < Am? + q^ 
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where ki and ^2 are determined by the condition ^- — " 
The intraband part of the susceptibihty thus reads 



> 1. 



g: 



Re 



(^)+-gn( 



167r 



sign 



2U3 



167r 



Gn 



+ / 2m- 



-G 




-G 

G^(^) 

-G^ 



G 



-G+ +G+ (?i^±^' 




lA 

IB 
2A 
2B 

3A 
3B 
4A 
4B 
5B 



Adding the interband part from above yields the final result given by Eq. (10 1. 



Appendix B: Relation between current and density correlation function 



We define the four-current J' 



The four-current correlator can then be written as 



9mX 



hA 



dte'"*-o*{(Mq,t),J"'(-g,0)])o 



Kg, 0, ^"(-9,0)]) J 



1 



Aujt—Ot 



dte 



itot — Ot 



QJ 



/5(q,t) + ig-i(g,0,(i + «)^"(-g,o) 



1 



{[p{q,0),J''{-q,0)])o.. 



where in the second line we used the continuity equation. This results in 

hA ^ 
1 

hA 



1 XjjM^^)i =^Xpj,iQ,^)q + 
1 



hA 



Pq,Q -J-q 
Pq,Q J-q 



The second term on the right-hand side vanishes while the third term needs special attention. For a translational 
invariant system, i.e., q'^Xj^j^ (q, a;)g' — q'^xfj^ one finally gets Eq. ([5|. In the ?^DEG, the last term is exactly canceled 



by the diamagnetic contribution, i.e., q^x^'^'^" — 0, ^-^d thus {d^J^)f^ = 0. In our Dirac model, gauge invariance is 
broken because of the cutoff in the valence band. This can be seen, for example, by lim^_»o ^^(g, 0) 7^ 0, which is 
unphysical, as a longitudinal static vector potential cannot induce a current. As stated in Ref.liS, taking into account 
the full Brillouin zone leads to the cancellation of the commutator by a diamagnetic contribution (which is absent in 
the linearized model) and thus to q^X^'^'^"qi' = 0- 
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